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Abstract: If, as is commonly accepted, the colour-singlet, 'preconfined', perturbative 
clusters are the primary units of hadronization, then the electric charge is necessarily com- 
pensated locally at the scale of the typical cluster mass. As a result, the minijet electric 
charge is suppressed at scales that are greater than the cluster mass. We hence argue, and 
demonstrate by means of Monte Carlo simulations using HERWIG, that the scale at which 
charge compensation is violated is close to the mass of the clusters involved in hadroniza- 
tion, and its measurement would provide a clue to resolving the nature of the dynamics. 
We repeat the calculation using PYTHIA and find that the numbers produced by the two 
generators are similar. The cluster mass distribution is sensitive to soft emission that is 
considered unresolved in the parton shower phase. We discuss how the description of the 
splitting of large clusters in terms of unresolved emission modifies the algorithm of HER- 
WIG, and relate the findings to the yet unknown underlying nonperturbative mechanism. 
In particular, we propose a form of as that follows from a power-enhanced beta func- 
tion, and discuss how this as that governs unresolved emission may be related to power 
corrections. Our findings are in agreement with experimental data. 
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1. Introduction 

Hadronization is one of the most poorly understood aspects of QCD. Our knowledge is 
currently limited to models ||l|, and estimations of the power-suppressed corrections Q 
to observables that have sensitivity to soft physics. 

Of the little that is known about the dynamics of hadronization, colour preconfinement 
1^, which is a general property of perturbative QCD, is often regarded as a plausible 
starting point. 

Colour preconfinement is a theorem, which follows from perturbative QCD, that states 
that in the course of the evolution between the hard scale Qhard and the cut-off scale Qq 
that results in the formation of a perturbative parton shower, the quarks and gluons become 
organized in colour-singlet 'clusters', whose mass is of order Qq and is independent of Qhard 
in the limit of large Qhard • 

It has been proposed |Q] that the clusters so produced participate independently in 
hadronization. If so, and provided that the physical cut-off scale to the parton shower is 
found to be small compared with Qhard, the nonperturbative contribution to jets should 
not dramatically disturb the properties of the perturbative parton shower. 

Models of hadronization based on colour preconfinement, notably the Monte Carlo 
event generators HERWIG ||l| and PYTHIA have been found to agree, to within cred- 
itable accuracy, with experimental data. When carrying out this comparison, since the 
dynamics of hadronization is little understood, it is necessary to tune some parameters 
that are related to hadronization. 
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One of the key parameters is the cut-off scale Qq, above which it is appropriate to 
apply the perturbative parton shower evolution and below which nonperturbative dynamics 
dominates. 

By means of experimental tunes it has been discovered that the cut-off scale is quite 
low, 0(1 GeV). Thus the typical cluster mass is also small. 

Despite the success of the models of hadronization based on preconfinement, the ex- 
istence of these preconfined clusters hitherto lacks direct experimental evidence^. In view 
of this, we turn our attention to the phenomenon of local charge compensation [^. We 
demonstrate that colour preconfinement naturally leads to local charge compensation. We 
advocate the measurement of minijet charge rather than distribution in rapidity space 
which has been traditionally considered. This change of observable allows us to relate the 
scale of charge compensation to the scale of colour preconfinement. 

In HERWIG, the perturbative clusters that remain large at the end of the parton 
shower phase are split according to a parametrization. We point out that because the 
mass distribution of the clusters is sensitive to soft emission that is considered unresolved 
in the parton shower phase, this cluster-splitting dynamics may be rephrased in terms of 
unresolved emission governed by a modified low energy running strong coupling. By com- 
paring the result with the default cluster mass distribution of HERWIG, we can estimate 
the energy scale involved in the cluster-splitting phase, or the shape of the modified as, 
and derive a possible phenomenological distinction. 

The cluster splitting energy scale thus established may be interpreted loosely as the 
scale of 'emission before confinement'. We discuss one explicit interpretation where en- 
hanced g ^ qq splitting modifies the running of as- The resulting form of low-energy as 
agrees with our findings with physically acceptable parameter values. We then discuss how 
this as relates to the part of the power corrections that is due to soft gluon emission. As 
the as has complex poles, there is ambiguity associated with its analytization. 

This paper is organized as follows. We first introduce the concept of colour preconfine- 
ment and the resulting local compensation of charge. We proceed to define the relevant 
observables. We carry out simulations using the HERWIG Monte Carlo event generator 
and compare the numbers against those of PYTHIA. We discuss the cluster-splitting proce- 
dure in HERWIG that affects this observable, and consider possible physics interpretations 
and consequences. Conclusions are stated at the end. 

2. Colour preconfinement and charge compensation 

In the original proposal of Amati and Veneziano [Q], units composed of perturbatively 
emitted quarks and gluons become colour-singlet clusters. Obviously the 'ends' of these 
clusters are defined by quarks. However, as the emission of quarks in a parton shower is a 
relatively rare event, the resultant clusters can be quite large even though the mass is still 
determined by Qq in the limit of large Qhard- In the Lund string model of PYTHIA, these 
quark-gluon systems are regarded as 'kinked strings' that subsequently decay into hadrons. 

^However, we note that the minijet structure of jets has been studied in the context of cluster hadroniza- 
tion in ref. M 
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An alternative, and more economical, approach adopted in the cluster hadronization 
model 1^ of HERWIG is to introduce forced, 'nonperturbative', splitting of gluons at the 
end of the parton shower, such that all clusters are quark-antiquark colour-singlet dipole 
systems. However, even in this case, there often remain in the end some clusters that 
are considered too large to be nonperturbative objects. These clusters are then split by a 
power distribution in terms of the masses of the decay products. We shall discuss this in 
more detail in sees. Q and ^. The small clusters at the end decay isotropically into hadrons 
according to the phase space weight. 

We may therefore say that the primary units of hadronization are, in the former picture, 
perturbative colour-singlet systems, whereas in the latter, all colour connected two-parton 
systems, or in other words, dipoles. 

In any case, if the physical cut-off scale is at <3o(physicai) j then for all artifical cut-off 
scales Qo greater than Qo(physicai) j emitted partons are organized into clusters with mass 
of order Qq. Therefore if we take an observable that is insensitive to physics below Qo, we 
would be measuring the properties of clusters with mass of order Qo- 

The electric charge of a cluster is always or ±1 because a cluster is in effect a quark- 
antiquark system. Following our reasoning above, the charge is or ±1 at all scales above 
Qo(physicai) ) such that if a cluster at a certain scale is composed of a number of smaller 
clusters at a smaller scale, the electric charges of these smaller clusters must mutually 
cancel. 

Hence electric charge is 'locally' compensated if colour is preconfined, rather than 
increasing in proportion to the square root of the number of charged tracks belonging to 
the cluster as would be expected if the charges were uncorrelated. It is worth noting that 
the converse is not necessarily true. In particular, in PYTHIA, although the strings may 
be quite large, their decay proceeds by the creation of quark-antiquark pairs, such that 
the scale at which charge is compensated is in general small compared to the mass of the 
strings. 

At scales lower than (5o(physicai) ) on the other hand, if nonperturbative effects are dom- 
inant, there is no reason to expect that charge is locally compensated. This is illustrated 
in fig. |l[ 



3. Charge compensation observables 

In practice, it is not possible to define an exclusive observable that is absolutely independent 
of physics below a certain arbitrary scale. However, we would physically expect that a 
minijet charge, if appropriately defined, can minimize the contamination. 

The sensitivity to physics below the scale defined by the minijet resolution variable 
depends on the algorithm that is used to combine tracks or objects consisting of tracks. 
We expect that based algorithms exhibit the most desirable properties since the 

resolution variable defined between two objects is, in principle, the energy scale that is 
required for a splitting that creates the two objects. In view of this, let us as our default 
procedure adopt the Cambridge |1C] algorithm. This algorithm adopts as the test variable 
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hard ^'^^ ^'^d 
process parton shower hadronization 
scale phase phase 

(charge not compensated) 

Figure 1: An illustration of the violation of local charge compensation during hadronization in 
the colour preconfined picture (HERWIG). Neighbouring charges do not necessarily cancel. 

between objects (ij) a variable that is essentially k'^: 

yij oc min(£'i, Ejfvij. (3.1) 

The normalization factor is taken to be l/E'^-^^ where -Bvis is the total visible energy. The 
quantity Vij, which also serves as the ordering variable, is defined by: 

Vij = 2{1- cos Oij). (3.2) 

The jet construction algorithm is iterative and the pair (ij) with the smallest value of Vij 
that satisfies yij < ycut is combined. In addition, for a pair which does not satisfy yij < ?/cut 
but have smaller values of Vij, the object with lower energy is 'frozen' by regarding this 
object as a jet. This last point marks the distinction between the Cambridge algorithm 
and the 'angular-ordered Durham' algorithm [p^ . 
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We now define the minijet charge as 



Qminijct(ycut) — ^ ^ Qtrack- (3-3) 

The summation is over the tracks that belong to the minijet under consideration. It is also 
instructive to define a charge that is weighted by powers of the track three-momenta, as: 

Q/ \ Qtrack (Ptrack ' Pminijet) A\ 

minijet (,2/cut ) /^j — ^^"7 yi j 

/ y vPtrack ' Pminijet j 



SO that comparison could be made with literature IQ. The k ^ limit of Qminijet(ycut5 ^) 
differs from the unweighted charge Qminijet(?/cut) by a factor of track multiplicity. 

Using the charge thus defined, we may further define the following quantities. First, 
the average minijet charge in an event is defined as: 

< |Qmimjct(ycut, ^)\ >= ^'Q'^°'j'-t(^-"^'^)l . (3.5) 

^minijets 

The measurement of the average minijet charge, measured and averaged over a sample of 
events, can confirm that charge is locally compensated. On the other hand, it would also 
be interesting to investigate to what extent the charge is limited to zero or ±1. We hence 
define the ratio: 

Rm>i = (3.6) 

Trminijets 

An alternative possibility would be to consider the average of some powers of minijet charge. 



4. Result of simulations 



We present the result of HERWIG simulations, and later provide a comparison with results 
obtained using PYTHIA, for the observables defined above. Unless stated otherwise, the 
simulations are at the hadron level and the default values in HERWIG 6.500 and PYTHIA 
6.215 are adopted for the parameters affecting hadronization. The HERWIG parameters 
are: 



RMASS(13) = 0.75, CLMAX = 3.35, CLPOW = 2.0, (4.1) 

RMASS(l) = RMASS(2) = 0.32, (4.2) 

PSPLT(l) = PSPLT(2) = 1.0, (4.3) 

CLDIR(l) = CLDIR(2) = 1, (4.4) 

CLSMR(l) = CLSMR(2) = 0. (4.5) 



Out of the parameters listed above, the most important are the following. The effective 
gluon mass RMASS(13) (in GeV) sets the cut-off scale for the parton shower, CLMAX (in GeV) 
sets the maximum allowed cluster mass, and PSPLT is the power for the mass distribution in 
the splitting of the clusters that have masses exceeding CLMAX. The cluster mass distribution 
before and after the cluster splitting is shown in fig. ^. For comparison, we also show the 
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Figure 2: The HERWIG cluster mass distribution. The curve 'HERWIG defauh' is generated by 
the default parameter settings. The curve 'large clusters' corresponds to the choice RMASS (13) =1.5 
GeV, CLMAX=5 GeV. We also show the result of having cluster splitting forbidden by setting an 
arbitrarily large CLMAX. The peak at about 5 GeV is due to clusters containing the bottom quark. 



distribution corresponding to RMASS(13)=1.5 GeV, CLMAX=5 GeV, which results in larger 
cluster masses than the default. 

To clarify, when the parton shower is terminated by the effective gluon mass, clusters 
are first formed by splitting the gluons into uu or dd. When they are heavier than ap- 
proximately CLMAX, clusters are split by generating an additional uu or dd pair. This is 
performed by uniformly generating mass raised to the power of PSPLT. At the end, clusters 
are decayed according to the phase space weight into the available hadrons. 

The other parameters, which are less important for our purpose, are the down- and up- 
quark masses RMASS(l), RMASS(2), the maximum cluster mass parameter CLPOW and the 
parameters CLDIR and CLSMR that define the extent to which the cluster decay remembers 
the direction of perturbatively produced quarks. 

For the main part of this section, we consider the simulation of the e~e~^ qq hard 
subprocess at the pole, ^/s = M^o = 91.188 GeV, such that kTcnt = Af^o -^ycut • Raising 
the energy raises the overall multiplicity and therefore raises contamination due to the mis- 
identification of charged tracks when constructing minijets. The initial/final state photon 
radiation is switched off in both HERWIG and PYTHIA. There is no detector simulation. 
The sample size is 10000 events except where otherwise stated. 

We first show the minijet multiplicity in fig. |^ Towards large j/cut we see that it tends 
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Minijet multiplicity versus resolution parameter 
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Figure 3: The minijet multiplicity plotted against the Cambridge algorithm resolution parameter 
j/cut- Sample size is 1000 events. The error due to finite statistics is small. 



to 2, whereas for small t/cut it begins to saturate as the number of tracks is finite. A ?/cut 
value of 10^^ corresponds to a kxcut of 0.91 GeV and hence of order of the cluster mass 
scale when using the default parameter values in HERWIG. 

In the same figure, we also show the result of adopting large clusters (RMASS(13)=1.5 
GeV, CLMAX=5 GeV). The difference between the two curves is slight, and the physics 
behind the difference is unclear. The minijet multiplicity itself is therefore not a sufficiently 
good observable for studying hadronization. 

We now turn our attention to the average minijet charge. We carry out our simulation 
again for the HERWIG default and modified parameter values. 

The result is shown in fig. ^ The Monte-Carlo sample size is 10000 events, but even 
with 0(1000) events, the statistical error is small enough to allow comparison of the two 
curves. We now find that in contrast to the minijet multiplicity, we have an observable 
whose behaviour can be directly related to the physics of hadronization. In particular, the 
peak position of the average minijet charge observable, where local charge compensation is 
maximally violated, is close to the cluster mass scale. Far below this scale, charge becomes 
that of individual hadrons, whereas above this scale, local charge compensation ensures 
that the minijet charge does not increase arbitrarily. 

Another point that is worth noting is that the minijet charge at large ycutj even in the 
two-jet limit, is fairly sensitive to the cluster mass. If the scale of hadronization is large, 
minijet charge is contaminated by the reshuffling of charged tracks amongst minijets. An 
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Average minijet charge versus resolution parameter 
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Figure 4: The mean minijet charge plotted against the Cambridge algorithm resolution parameter 

2/cut- 



alternative way of seeing this is through plotting the fraction of minijets that have charge 
exceeding 1, as shown in fig. |5[ 

The relation between the minijet charge and the multiple-charge fraction can be un- 
derstood as follows. Let us first consider the charge of a quark or a gluon jet in the ideal 
case where all clusters emitted from the originating parton is assigned to this jet. In this 
case, the net charge of this jet has to be equal to or ±1 with equal probability. This is 
because the charge is determined by the typically nonperturbative splitting g ^ qq that 
is closest to the hard process. All other emissions conserve the jet charge. This produced 
quark has equal probability of being a u or a d. For the case of a quark jet, for instance, 
regardless of whether the originating parton is up-type or down-type, the net charge of the 
jet is therefore always or ±1 with equal probability. Hence the 'perturbative' jet charge 
should average to 0.5 so long as the charges of different jets are uncorrelated, or perfectly 
correlated in the case of the two-jet limit. 

Now let us consider the contamination of this jet with at most one singly charged 
track, whose charge is uncorrelated with the charge of the jet, with equal probability p for 
each charge, +1 or —1. We can show that the probabilities for the total jet charge are now 
modified to: 

1 - 2p » 1 



IQI=o 



-2p p 
2 2 



p 



Pi 



IQI=i 



2p 1 - 

— H 

2 2 



2p 



(4.6) 
(4.7) 
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Multiple charge ratio versus resolution parameter 
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Figure 5: The fraction of minijets that have charge exceeding 1, plotted against the Cambridge 
algorithm resolution parameter ?/cut- 



(4. 



_ P 
1QI=2 - 2' 

Summing the charge times probability for each case, we obtain the expectation value for 
the net charge as p+1/2. On the other hand, the multiple-charge fraction is given by -F|q|=2 
above, such that we obtain the relation between the minijet charge and the multiple-charge 
fraction as: ^ 

< IQminijctl >= 2 + 2-^IQI>i- (4-9) 
Comparing figs. ^ and ^, we see that this relation is satisfied very well in the two-jet 
limit. For lower values of ycut) as the nonperturbative effects become more dominant, the 
relation is increasingly violated, although the general behaviour is still in accord with the 
expectation from eqn. (|4.9| ). 

The assumption here of contamination due to at most one charged track is sufficient 
for our discussion, but for the sake of completeness, let us derive in outline the case with- 
out restriction on the number of uncorrelated contaminating tracks. We first define the 
generating function for the jet charge by: 



1 1 u 

4^ 2 4 



exp 



P 



- -2 + u 
u 



(4.10) 



We may obtain this exponential form as the limit of binomial distribution due to infinitely 
many contaminating tracks, which is a direct generalization of eqns. ( [4.6D -( |4l^ ). We recover 
eqns. ( |4.6D -( p^ in the limit of small p. 
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Bessel functions of the first kind are defined by the generating function: 



exp 



z ( 1 



J2 *"'^«(^)- (4-11) 



After a few elementary algebraic manipulations, we can equate the two generating functions 
to obtain: 



\i^-^Jn-i {-2ip) + U'^Jn {-2ip) + J„+i {-2ip) 



PQ=n = e-2P 

Bessel functions can be expanded in a series: 



(4.12) 



~ 2«+|m| 

{-2^p) = ^ J (4.13) 



such that corrections to eqns. ([4.(]|)-( p^ ) can be obtained systematically as an expansion 
in p. Furthermore, we observe that by successively operating on eqn. ( 4.1C| ) by u{d/du), we 



may obtain the expectation value for charge raised to any even integer power. For instance, 
operating twice with u{d/du) we obtain: 

< Qminijet >=\ + ^P- (4.14) 

We can estimate the contamination probability p naively as the ratio of the cluster 
mass against the jet energy multiplied by the cluster multiplicity: 

p~^x#d.,. (4.15) 

The cluster multiplicity can be estimated from fig. ^. As the cluster size increases, the clus- 
ter multiplicity decreases and so there is some cancellation between the two contributions. 
This naive estimate can be compared with the two-jet limit of fig. |5| and we see that both 
the cluster mass dependence and the overall magnitude is described well. Thus the net jet 
charge in the two-jet limit by itself already provides a measure of the cluster mass scale. 

Having said this, it is important to test to what extent the contamination is a feature 
of nonperturbative dynamics and not a defect of the jet algorithm, i.e., misclustering. To 
this effect, we repeat the analysis using different jet algorithms and show the result in fig. |6|. 

The JADE algorithm |ll, that uses the invariant mass rather than kT as the resolution 



variable, has peak at higher y^ut than the /cj^-based algorithms. One reason is the different 
kinematics, namely that the kx between two tracks is always smaller than their invariant 
mass. Another reason is the increased mis-clustering at high j/cut where the dynamics is 
perturbative, and reduced mis-clustering at low ycut where the dynamics is nonperturbative. 

Before discussing this point in more detail, we turn our attention to the three remain- 
ing, /cT-based, algorithms. The mis-clustering at moderate ycut is largest for the Durham 
algorithm and is slightly better for the angular-ordered Durham algorithm that was men- 
tioned in sec. ^ At very low ycutj where angular-ordering is no longer a feature of the 
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Average minijet charge versus resolution parameter 
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Figure 6: A comparison of tlie minijet charge obtained using different jet algorithms. 



relevant dynamics, the Durham algorithm actually fares better than the angular-ordered 
algorithm. 

At large j/cut the minijet charge due to the four algorithms converge. In particular, 
there is little distinction between the large ?/cut behaviour of the Cambridge algorithm and 
the angular-ordered Durham algorithm from which the former algorithm is derived. Hence 
it seems reasonable to suppose that the minijet charge measured using the Cambridge algo- 
rithm represents the limit to which mis-clustering could be suppressed in the perturbative 
large ycut region. In this sense, the ansatz made earlier that the contamination probability 
p has a mainly nonperturbative origin is reasonable. 

On the other hand, towards very low i/cut) as we have observed above, the invariant 
mass, as in the JADE algorithm, may become a better resolution variable than kx- Hence 
the measurement of minijet charge using these two algorithms could provide complementary 
information. 

This suggests the introduction of a new class of scale-dependent jet algorithms, which 
combines the low yij behaviour of JADE with the large yij behaviour of /^T-based algo- 
rithms. As an example, we propose a 'preclustered' scheme, where tracks are first combined 
using the JADE algorithm up to some point given by the resolution parameter ?/had) after 
which the tracks are combined using a /c-r-based algorithm with resolution parameter 2/cut- 
The case that uses the angular-ordered Durham algorithm is shown in fig. |^. The preclus- 
tering cut-off yhad is the normalization times (0.5 GeV)^. This value is the most successful 
that we have been able to find when using HERWIG with default parameter values. 
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Figure 7: Minijet charge in the 0.5 GeV 'preclustered' angular-ordered Durham algorithm com- 
pared with the JADE and angular-ordered Durham algorithms. 



There is marked improvement, in terms of the peak height and large ycut behaviour, 
compared with either JADE or the angular-ordered Durham algorithm. There is also 
improvement compared with the Cambridge algorithm. 

We found that there is less marked improvement when preclustering is followed by 
either the non-angular-ordered Durham algorithm or the Cambridge algorithm, and there 
is no visible improvement when the JADE preclustering phase is also angular-ordered. 

The preclustered curve remains large at low ycut compared with the JADE and angular- 
ordered Durham curves, and in the low ycut hmit tends to the value obtained using the 
JADE algorithm at the point ycut = yhad- 

The sensitivity of the minijet charge to the parameters affecting hadronization, as 
well as to the details of the jet algorithm, suggests that because of our limited knowledge 
about the dynamics of hadronization, the prediction of any existing Monte-Carlo event 
generator can not be completely trusted when charge distribution is concerned. On the 
other hand, the violation of local charge compensation during the hadronization phase 
is presumably a universal phenomenon in all models of hadronization that are based on 
preconfinement. Thus the general behaviour of the minijet charge observable, that peaks 
around the hadronization region, is a concrete prediction for this class of models. 

To demonstrate this point, in fig. we show the numbers obtained using PYTHIA. The 
overall behaviour is fairly similar, except the behaviour at very small ycut; and there is an 
unwelcome dependence on whether the decay is matrix-element corrected, although we 
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Figure 8: A comparison of HERWIG and PYTHIA for the mean minijet charge. The Unes PYTHIA 
(0) and PYTHIA (2) correspond to the uncorrected option and the second-order matrix-element 
corrected option respectively for the Z'^ decay in PYTHIA. 



have found that the difference between the first-order and second-order corrected options 
is small. We have confirmed that the corresponding difference between the matrix-element 
corrected and uncorrected options is almost negligible in HERWIG. The effect of adding 
the matrix-element correction in HERWIG is merely to raise the charge slightly at larger 
values of ycut as expected. 

As mentioned before, although the PYTHIA picture is that of string fragmentation 
and the mass of the string may be large to start off with, the scale at which charge is 
compensated is much smaller because the fragmentation proceeds by the creation of quark- 
antiquark pairs. We see from the figure that the peak position is roughly the same as that 
of HERWIG. 

In fig. m we plot the momentum squared weighted minijet charge, which corresponds 
to K = 2 in eqns. (0) and (3^). The bottom quark contribution is plotted on a separate 
curve in order to allow comparison with ref. [12| in the two-jet (large ycut) limit. 

The general behaviour is quite different from the unweighted charge of fig. ^. In the 
low ycut region, the weighted charge tends to the same value as the unweighted charge. 
This is interpreted as the average hadron charge. In the high ycut region, the charge is 
much lower compared with the unweighted charge. This is because the observable suffers 
less from 'contamination' due to soft tracks. 

We have also found that for a large and increasing k (10 ~ 50) and in the two-jet limit. 
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Figure 9: The momentum-squared weighted minijet charge, plotted against the Cambridge al- 
gorithm resolution parameter j/cut- In addition to the HERWIG default (solid line) and the large 
cluster (dashed line) settings, we also show the bottom quark contribution by the dotted line. 



the charge slowly increases in the region 0.5 ~ 0.6. In the limit of large k, the momentum- 
weighted charge tends to the charge of the highest energy hadron in each jet. From figs. ^ 
and P we estimate the average hadron charge to be slightly above 0.6, in agreement with 
the above finding. 

For moderate values of k, it is not reasonable to assume that only the highest energy 
hadron contributes. In this case, the charge is estimated as follows. Let us assume that 
the highest energy hadron originates from the isotropic n-body decay of a heavier object. 
If this object is a cluster, n is on average ~ 4, whereas if it is a heavy intermediate hadron, 
n is typically 2. It is a reasonable approximation to assume that this object has charge 
or ±1, and it is also reasonable to assume that the charge of each of the decay products 
is also or ±1. Then the expectation value for the momentum- weighted charge of this 
object, when it decays into #tot objects out of which #ch are charged, is simply: 

We assumed massless kinematics. For small #tot or #ch5 this multi-dimensional integral can 
be evaluated analytically, but the general case is presumably best evaluated numerically. 

Let us now adopt the viewpoint that the momentum-weighted charge is mainly due 
to the heavy intermediate hadron which decays by 2-body kinematics to give the highest 
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energy hadron in the jet. The other contributions are expected to mainly affect the de- 
nominator of eqn. ( 4.16| ) by increasing it, as uncorrelated contributions to the numerator 
cancel on average so long as they are not too large. 

For the ^^tot = 2 case and k = 2, we obtain 0, 0.5 and log 2 = 0.693 respectively 
for = 0,1,2. Disregarding the contribution from other hadrons, we can estimate the 
momentum-weighted charge as some weighted average of these three numbers. We note 
that out of the three numbers, only the #tot = 1 case corresponds to nonzero overall charge. 
Hence one possible estimate would be 0/4 + 0.5/2 + 0.693/4 = 0.42. This is somewhat large 
compared with the large ?/cut end of fig. ^, such that we see that the contribution from other 
hadrons is not negligible but is under control. 

Our numbers for the mean momentum-weighted jet charge can be compared with the 
hemisphere charge separation measurement in ref. [12|. In the two-jet limit, the momentum 
component along the jet axis would not be very different from the momentum component 
along the thrust axis, such that we expect that the mean weighted charge presented here 
would correspond to a half of the charge separation of ref. However, a quick compari- 
son, for example with their fig. 2, shows that there is nearly a factor of two difference. The 
hemisphere charge separation is much smaller than is expected from HERWIG. This may 
be due partly to the 'detector effects' included therein. Furthermore, the numbers due to 
JETSET presented in their tab. 4 has much greater flavour dependence than in our work, 
even in the k ^ oo limit. Further study is desired in order to elucidate the nature of the 
differences. 

For comparison, in fig. IC, we show the result due to PYTHIA (JETSET is now part 
of PYTHIA). In the two-jet limit, the two generators give roughly equal numbers, such 
that the difference between HERWIG and ref. [|l^ is most likely not due to the difference 
in the generator. For small Ucut, the difference that is already visible in fig. ^ is magnified 
and the PYTHIA numbers show a peak at a few times 10~^. 

The numbers shown so far have been obtained for the simulation of e^e~ collisions at 
the pole. On the other hand, considering the high energy experiments that are currently 
in plan, it is interesting to consider how our procedure and results are modified for collisions 
in a hadron collider environment, where there are contaminations from emission near the 
beam axis from the initial-state parton and also from the soft underlying event. 

We generate QCD 2 — > 2 hard scattering events in pp collision at 2 TeV in the central 
region. We require for both outgoing partons rapidity \y\ < 0.5 and moderate transverse 
momentum 30 GeV < pt < 50 GeV. Events are generated with and without soft underlying 
events simulated using the HERWIG default procedure and with default parameters. 

Minijets are defined using the /c-r-type jet algorithm of KTCLUS |13] in this case as 
follows. We first define macrojets by selecting a large enough 7/cut value such that exactly 
two jets remain uncombined with the beam. We then construct subjets constituting these 
two macrojets at new values of i/cut- These subjets are regarded as the relevant minijets. 
In effect, we are studying the charge structure of tracks within the two macrojets. The 
normalization factor in eqn. (|3.1| ) for ?/cut is taken to be 1/M^o rather than l/E"^^^, such 
that the scale matches that of the e^e~ simulations. 

The result is shown in fig. As the jet algorithm of KTCLUS is not angular-ordered. 
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Momentum squared weighted charge versus resolution parameter 
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Figure 10: A comparison of HERWIG and PYTHIA for the momentum-squared weighted minijet 
charge. The uncorrected decay option is adopted for PYTHIA. 



we should compare the result against the Durham algorithm shown in fig. When the 
soft underlying events are suppressed, the low ycut behaviour is similar to that of fig. |6|, 
but the peak has weakened and the minijet charge remains large at large ycut- From this 
curve alone, it seems that the best that can be achieved in a hadron collider environment 
is a fit with Monte Carlo simulations. On the other hand, further refinement of the track 
selection procedure may improve the prospect. The additional consideration of the JADE- 
type algorithm is another possibility. 

The effect of adding soft underlying events is quite marked. However, we note that the 
default procedure for soft underlying events in HERWIG violates local charge compensation 
such that the numbers may be regarded as a pessimistic estimate. In the alternative 
HERWlG-based simulations of ref. |14], as well as in PYTHIA, the underlying events are 
generated as multiple parton scattering and so charge is locally compensated per each 
scattering. Having said this, whether charge is locally compensated in soft underlying 
events at hadron colliders is an open question which merits further investigation. 



5. Interpretation in terms of cluster dynamics 

As mentioned in the previous sections, the mass of the perturbative colour-preconfined 
clusters is of order of the cut-off scale. Put in other words, emissions which are normally 
considered unresolvable in the parton shower in fact affects the mass of clusters. If so, it 
is perhaps natural to consider the splitting of large perturbative clusters, that occurs in 
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Average minijet charge versus resolution parameter 
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Figure 11: Minijet charge in pp collision at 2 TeV. QCD 2 — > 2 scattering events are generated in 
the region |y| < 0.5 and 30 GeV < pt < 50 GeV, with and without soft underlying events. Only 
tracks belonging to the two macrojets are considered. ?/cut is normalized by l/M|o. 



HERWIG, in terms of an extension of perturbative dynamics using a modified running as, 
which describes emission below the cut-off scale 'before confinement'. 

There are several arguments in favour of this modification. First, in the default HER- 
WIG picture, the large clusters are split by a string-like mechanism, by creating quark- 
antiquark pairs in the string vacuum. On the other hand, the momentum direction of the 
quark and the antiquark constituting the decaying cluster is conserved as if the cluster is 
still an object composed of a perturbative quark and antiquark, unlike in the subsequent 
cluster decay which is isotropic. Second, and in principle, this substitution could lift the 
scale dependence in HERWIG between the parton shower phase and the cluster splitting 
phase such that the ambiguity due to the incomplete description of hadronization is in prin- 
ciple reduced to still lower scales. Third, the necessity to invoke an uncalculable string-like 
mechanism is removed, and we would instead have a description in terms of a modified 05, 
whose properties could be subjected to further study. 

Before going into a discussion of the meaning of this as, let us first consider the Sudakov 
form factor due to the left-over 'unresolved' emissions. The addition of an 'unresolved' 



gluon to an angular-ordered parton shower is illustrated in fig. 12, where we consider the 
emission of an unresolved gluon in between rapidities yo and 7/1, where yo and yi are the 
rapidities of the gluons which have been emitted in the course of the resolved parton shower. 
The emission probability of a gluon in a phase space element {dQ'^, dy), where and 
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Figure 12: The addition of an 'unresolved' gluon to an angular-ordered parton shower in between 
colour connected partons (gluons). y is the rapidity of the emitted parton. 



y are the and rapidity of the emitted gluon respectively, in the soft limit, is: 



"■-'emission — Uy 

vr 



The colour factor C is Cp for emission from a quark and Ca for emission from a gluon. As 
we are considering cluster dynamics and the gluon is effectively a quark-antiquark object 
in this respect, the colour factor is C_f, or equally well Ca/2, at leading order in 1/Nc- 

The equivalence of this expression with the ordinary angular-ordered parton shower 
1 16] in the soft limit is demonstrated as follows. Let us consider the emission probability 
from a quark in the splitting Qa — > Qbgc- 

dPemission{qa ^ Qbgc) = ^^0!siQ^)Pqq{zb). (5.2) 

t 2tt 

Here i = E'^{1 — cos Obc)- Ea is constant until an emission takes place, such that di/i = 
dffi jd'^ = 2dy in the soft and the usual small angle limit. In the soft limit, the splitting 
function Pgq{zb) = 2Cf/zc, such that we have: 

dPemissioniqa ^ qbOc) = 2dy'^as{Q'^)^^ ■ (5.3) 

zvr Zc 

Now replacing dzb by dzc and using = 2iz1 in the soft limit, we see that the two 
expressions are equivalent. 

For an angular-ordered parton shower, the emission is ordered in y. The Sudakov form 
factor is: 

dQ^ , as{Q^)CF' 



Ag{y) = exp 



-dy- 



(5.4) 



Q2 TT 

We are only interested in emission below the cut-off scale, such that the integration 
factorizes out: 

Ag{y) = [exp{y - yo)]-'\ (5.5) 

where 

w — ^ — • ^'-'^ 
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The integral is finite if as tends to zero as a power as tends to zero. A possible physical 
argument would be to say that if colour is confined in the large distance limit, then any 
as with some physical interpretation as the coupling of coloured objects can be defined to 
vanish as a power. 

We may thence carry on to estimating Iq in the above using models of low energy as- 
As an example, and not necessarily an appropriate one, we may consider the effective as 
of ref. [^5|, which is derived from the general consideration of power-suppressed corrections 
to event-shape observables. In this approach, the renormalon chain which gives rise to the 
power-suppressed corrections is summed using a dispersion relation. The effective as is 
defined to be analytic everywhere except along the negative real axis, and its infrared 
behaviour is connected to the power corrections through the analyticity of the observable 
under attention, or more precisely its characteristic function, as a function of small squared 
'gluon mass'. 

One difficulty when using this as to estimate Iq is that it remains finite as — > 
as shown in fig. |l^, such that its integral is not finite. In the dispersive resummation of 
the renormalon chain, this is not a problem as only the moments of as are related to the 
power-suppressed corrections. For now we content ourselves by cutting off at Q ~ 0.1 GeV 
and integrating up to the HERWIG cut-off at around RMASS(13)=0.75 GeV. We obtain 
Jo ~ 1- 

We shall discuss the above points regarding the behaviour of low-energy as in more 
detail in sec. ^. 

It is a simple matter to relate the emission dynamics to the kinematics of cluster 
splitting. Let us consider the splitting of a dipole with mass Mq into two smaller dipoles 
with masses Mi and M2. The kinematics has already been calculated in the context of the 



dipole cascade formalism in ref. |17]. As the integral is Lorentz invariant, let us transform 
to the cluster frame. Neglecting quark masses, rapidity in this case is: 

y = log(M2/Mi), (5.7) 
\y\ < log(Afo/Qcut). (5.8) 

Combining these with eqn. (^]^) and /q ~ 1, we may estimate the probability of there 
being no splitting. If the cluster mass is ten times greater than the cut-off, for example, 
the probability of the cluster surviving without emission is only 1%. 

In general, as can be inferred from eqns. (|5.5| ) and (|5.8| ), Iq governs the average number 
of smaller clusters that a large cluster splits into, viz.: 

<#.piit >~ /oAy « Jo log -2^. (5.9) 

^cut 

From fig. |2|, we can estimate the value of <#spiit> that leads to the correct cluster multi- 
plicity. We estimate this number to be ~ 3 and Mq to be ~ 10 GeV. This gives Iq ~ 0.5 
and thus /q ~ 1 obtained above is too large. 

= Pj^ on the other hand controls the typical cluster mass, and can be expressed as: 
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Effective low energy ttg versus energy scale 
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Figure 13: The effective as of ref. ||l5|, derived from the general consideration of power-suppressed 
corrections to event-shape observables. The effective as coincides with the three-flavour as in the 
high-energy limit. The three-flavour Aqcd is taken to be 0.25 GeV. We also show a Gaussian 
approximation, with arbitrary normalization. 



If we omit the consideration of the nonperturbative g ^ qq splitting that converts the 
dipoles into clusters, the physics of which we do not yet understand and which we will 
discuss in sec. ^, the cluster splitting is henceforth governed by the above equations. One 
should add a technical detail, that in HERWIG, the cluster mass that is considered has 
the constituent quark masses subtracted. 

We note that the default cluster splitting procedure of HERWIG violates the above 
expression, eqn. ( 5.10 ), for Q^, when the product of the masses of the two decay products 
divided by the decaying cluster mass is greater than the parton shower cut-off. 

Based on the above considerations, we carry on to creating a simple code which modifies 
HERWIG to perform cluster decay by the modified-as approach. However, there are two 
subtleties that need to be dealt with. 

First, if there is a cascade splitting of clusters, one must be careful not to double-count 
soft emission. The simplest method for avoiding double-counting would be to introduce a 
flag. On the other hand, when making a HERWIG implementation, as there is a cut-off 
on the maximum allowed cluster mass, this problem is not severe. 

Let us consider again the cluster mass distribution shown in fig. ^. The parton shower 
cut-off, which also signifies the scale of the low-energy cluster-splitting dynamics, is of 
order RMASS(13)=0.75 GeV, whereas the minimum cluster mass that is split is around 
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CLMAX=3.35 GeV. We see that in most of the cases, provided that eqn. (|5.1[l|) is satisfied, 
at least one of the decay products is fighter tfian tfie cut-off. Hence there would be no 
further emission from this lighter cluster, and by adopting an ordering that starts from the 
largest value of \y\, the problem would be resolved. 

Second, and related to the first point, the ordering in y is not strictly necessary. So 
long as the cluster mass can be considered large compared with the parton shower cut- 
off, the order in which the 'unresolved' gluons are emitted, in principle, does not matter. 
However, in a practical HERWIG implementation, because of the cut-off on the maximum 
allowed cluster mass and the incomplete description of physics below this scale, there does 
remain some dependence on the ordering. 

In particular, when using the above effective as, because /q ~ 1 is large, we have found 
that if emission is in the order of decreasing \y\ as suggested above, too many clusters are 
formed with small mass. This is unwelcome in the sense that these small clusters can not 
be considered perturbative objects. On the other hand, when /q ~ 0.5 as is preferred from 
multiplicity considerations, we have found that there is not much dependence on whether 
1 7/1 is ordered. In this case there is not a clear-cut connection between the presence of 
this ordering and, for example, the cluster mass distribution, such that for now, we may 
abondon the ordering in y and generate it uniformly. 

After the above considerations, we implement the following procedure, which is not 
precisely the procedure dictated by the modified-as approach, but more close to the HER- 
WIG default cluster-splitting algorithm: 

1. Terminate the parton shower by means of an artificial gluon mass RMASS(13)=0.75 
GeV, as according to the HERWIG default prescription. 

2. For clusters whose mass is greater than about CLMAX, or more precisely for a cluster 
with mass M made from quarks i and j: 

^CLPDw > cLMAX^L™ + (RMASS(i) + RMASS(i))^^™, (5.11) 
split the cluster, as follows. 



Rapidity is first generated, not by following eqn. ( p. 51) , but as a fiat distribution in 
the range suggested by eqn. ( |5.8D . This is correct in the limit where there is one and 
only one 'unresolved' emission. It is also correct in the case where emission is ordered 
in Q^, as in ref. |17|. In this case, the distribution of described below should be 



interpreted as a Sudakov form factor rather than as- 

is generated in between the parton shower cut-off scale and the lower cut-off scale 
for as, according to the distribution: 

The process is repeated until there remains no more cluster that satisfies eqn. ( 5.11J ). 



oc ^as{Q'). (5.12) 



Thus all clusters above the cluster-mass cut-off are split. The error due to disallow- 
ing some heavy clusters is negligible so long as the probability of splitting, or the 
parameter Iq, is sufficiently large. 
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HER WIG cluster mass distribution 




cluster mass /GeV 

Figure 14: The cluster mass distribution generated using the HERWIG default procedure, the as 
of ref. pa and a Gaussian approximation with Qmin = 0.1 GeV. 



In the above implementation, the normalization of as, which should affect multiplic- 
ity, is nevertheless irrelevant so long as we abandon y ordering. Hence one may as well 
substitute for low energy as by a Gaussian distribution. A useful substitution may be: 



oc Q exp 



{Q-Q 



^2 



mean ) 



2(AQ)2 



(5.13) 



such that Q is generated with mean Qmean and standard deviation AQ. Although this is 
a fairly economical parametrization, we may further reduce the number of parameters by 
setting the upper limit equal to the parton-shower cut-off. Thus we set: 

_ RMASS(13) + Q^in 

^jmean — 2 ' V^'^^) 

AQ = ^^ASS(13)-Q...^ 



The Gaussian distribution corresponding to Qmin = 0.1 GeV is shown in fig. 13. 

The cluster mass distribution due to this choice is shown in fig. |l^, where we compare it 
against the HERWIG default distribution and against the result of adopting the procedure 
above with the as of ref. [^]. When using this as, is generated in between 0.1 GeV 
and 0.75 GeV. 

We see that with Qmin = 0.1 GeV, although the distribution is close to the HERWIG 
default, the distribution is too skewed to the low mass end, and multiplicity is not large 
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Average minijet charge versus resolution parameter 
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Figure 15: The mean minijet charge plotted against the Cambridge algorithm resolution parameter 
j/cut- The sample size is 50000 events for each curve. 

enough. These two problems are mutually related, but presumably the more serious is the 
skew to the low-mass end, as the multiplicity is in principle affected by the normalization 
of as as noted above. 

Both problems can be rectified by raising Qrain^ and we have found that the correct 
enhancement of the high mass end is obtained at Qmin ~ 0.4 GeV. The individual identified 
particle yields are also affected by this cluster mass distribution. As the shape of the 
distribution is not identical with the HERWIG default, the yields would be slightly different, 
but as there is no reason to believe that the Gaussian distribution resembles reality, there 
would be no reason to believe that the resulting yields would be closer or farther to reality 
than the default HERWIG numbers. 

In fig. we show the average minijet charge. The numbers are stable with respect 
to small perturbation of the parameter Qmin- The result shows that the peak position in 
ycut has shifted slightly to the left compared to the HERWIG default numbers, and the 
peak height is smaller. As the cluster mass distribution itself is more or less unchanged, 
this modification should be interpreted as the result of having less probability of a large 
cluster splitting into two adjacent large clusters. Thus a modification of HERWIG that 
excludes cluster splitting that can be interpreted as due to high pT emission leads to the 
suppression of minijet charge. 

This argument also helps to understand why the PYTHIA minijet charge, without the 
matrix element correction, is larger than the HERWIG result, as seen in fig. M. 
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It is interesting to see whether event shape observables are also affected. We have 
computed sphericity and found that the modification due to the implementation of the 
new algorithm is slight. 



6. Interpretation of low energy as 

What we have learnt from the above simulation in terms of the dynamics of hadronization 
is that the scale relevant to the splitting of clusters is between about 0.4 GeV and 0.75 
GeV. In terms of the perturbative as, between about 0.5 and 1. 

What remains is to understand the origin of this structure. If we reduce the scale 
of splitting, the cluster mass distribution is too skewed towards the low mass end and 
multiplicity becomes insufficient. It would not be sufficient to propose, for instance, that 
only the highest pT emission affects the splitting of clusters. 

Let us consider a toy model, in which perturbative emissions have Sudakov form factor 
A((5^) and the confinement interaction has Sudakov form factor Ac{Q'^)- Then the prob- 
ability of emission 'before confinement' is given by exponentiating the emission probability 
times the survival probability and hence the following modified Sudakov form factor: 



Aefrective(Q^) = exp 



-Ac(Q 



exp 



g2 ^ 

—^dyAciQ ) 



(6.1) 
(6.2) 



Although the confining potential is not obtained at any finite order of perturbation theory, 
for the sake of discussion we may write confinement as an 0(a|) 'gluon exchange' process. 
Thus Ac(Q^) = 1 - 0{al) and hence we have a plausible explanation for the suppression 
of emission at lower scales where as{Q'^) becomes large. 

In terms of probability, this is equivalent to saying that although the emission of gluons 
at low has large probability, in the time that it takes to have this emission, there is 
greater probability that the cluster will be confined. Hence emission is suppressed. 

In terms of effective as, the effective as that controls 'unresolved' emission that splits 
the cluster is expanded in terms of perturbative as as as — 0{a^). However, the coefficient 
may depend on the cluster mass and the nature of the observable and hence this as is not 
necessarily universal. 

In the above model, we have assumed that the onset of confinement is sudden. This 
is also the case in HERWIG, where the 'confined', 'nonperturbative' clusters thus formed 
then decay isotropically. In reality, the soft gluon exchange interaction between the quark 
and the antiquark in a cluster is expected to enter earlier. Hence clusters would be 'rotated' 



as illustrated in fig. 16. We expect that this leads to the broadening of the peak in the 
minijet charge observable, although a quantitative prediction is beyond the scope of this 
study. 

Let us illustrate and concretize the above discussion of 'emission before confinement' 
by considering the dynamics of the g ^ qq splitting intrinsic to HERWIG which is not 
understood. If this splitting is physical, for instance if confinement enhances this splitting. 
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Figure 16: Perturbative and nonperturbative dynamics affecting perturbative clusters. 



it would be possible to define the corresponding beta function. Solving the renormalization 
group equation, we may define another as that may be regular in the infrared. This 
provides another viewpoint that is compatible with the above toy model in the sense that 
we may loosely identify 'confinement' with the g ^ qq splitting. 

We may define the effective strength for gluon splitting, that could be used to replace 
HERWIG's non-perturbative gluon-splitting, as follows: 

IICa 



(6.3) 



where /? = da^^ /dlogQ^ . Teg = n/Tfj in perturbation theory at the lowest order, but it 
becomes large at low scales for any parametrization of as that turns over. If we require 
that as vanishes as a power of such that its integral with respect to log is finite, 

we see that Tcs must diverge towards —>^ with the same power, i.e., asTcs tends 
to a constant. Although this power behaviour can not be obtained at any finite order in 
perturbation theory, it is precisely what one would expect as the form of a nonperturbative 
higher-twist contribution. We may write: 



/3(Q')=^pert.(«5)-^. 



Dropping the higher order perturbative contributions, we obtain as as: 

1 



asiQ' 



(6.4) 



(6.5) 



/3olog(Q7A2) + 2^/(nQ")' 

Here j3o = (33 — 2n/)/127r. Although the forced g ^ qq splitting in HERWIG only involves 
two fiavours, for the sake of consistency with the previous discussions, let us adopt n/ = 3 
for now. This as has the desired properties of integrability under log and remaining well- 
defined on the positive real axis, provided that A > /3oA"e^^ such that the denominator is 
positive definite. On the other hand, when either the denominator is positive definite or 
n > 4, the expression has complex poles. 
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If we require that this as is universal, we should require that n is large such that 
we do not introduce power corrections that are not predicted by the OPE formalism to 
quantities that are proportional to agiQ"^) in leading order such as the e"'"e~ total cross 
section. OPE predicts a term proportional to the dimension-4 gluon condensate, which 
is therefore oc l/Q^ by dimensional analysis, whereas the above as gives rise to terms 
(X [a^^)"^ /Q"^ for large Q, where a^^ is the perturbative as, such that the n = 2 case is 
forbidden. 

In the discussion of ref . [ 18 1 , even additional contributions proportional to 1/Q^ is dis- 
favoured, as this gives an additional 'nonperturbative' contribution which is not related to 
the region of small momentum flow, i.e. the region which is responsible for the condensates 
that give rise to the power corrections, in the corresponding Feynman diagrams. 

As an example, let us adopt n = 8 here, which is the same as the large behaviour 
of the effective as of ref. [^]. We set A = (1.8A)"'/3o and display the behaviour of this 
perticular choice in fig. 17. For reasonable choices of the parameter A, the low energy 
cut-off is large, in quantitative agreement with our earlier finding that Qmin ~ 0.4. On the 
other hand, the low behaviour of the aforementioned as is matched better for small n. 
In fig. 17, we also plot the n = 2 case for comparison, with A = (1.4A)"/3o- 

As stated above, the n = 2 case is forbidden and a 1/Q^ correction is disfavoured by 
OPE, but we claim that this argument does not apply to our 05 at n = 4, since what we 
are considering is a possibly universal correction to as arising from OPE, which therefore 
arises from the region of small momentum flow. 
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Figure 18: An illustration of the splitting of a 'kinked string' by the parton splitting g qq. 
Classically, the string 'pulling apart' the quark and antiquark approximates to a string with negative 
tension 'pushing apart' the partons. 



If the g ^ qq splitting is indeed enhanced by confinement, an intuitively appealing 
picture would be that the quark and the antiquark are 'pulled apart' by the large-distance 
string-like confining potential, as shown in fig. |l^. This would classically approximate to 
a string with negative tension 'pushing apart' the quark and the antiquark, and hence 
by dimensional analysis the contribution to the effective strength for gluon splitting is 
expected to be: 



cS 



UfTR 



(6.6) 



njTji = 1 or 1.5 depending on whether the splitting g ^ ss is allowed. Hence A ~ cr^. Let 
us test the validity of this parameter choice. First, the peak position in this as is given by 



the point at which the beta function vanishes in eqn. (6.4), i.e.: 



A — PoQ peak- 



Now from eqn. KW) and g^in = 0.4 GeV, Qpeak = 0.575 GeV. Thus A 



peak 



(6.7) 
0.237 



GeV^ is in good agreement with what one might expect for the string tension a ~ 0.2 GeV^ 



1 19]. The corresponding as is shown in fig. 17. 

The area underneath as and hence Iq in eqn. ( ^.5]) can be calculated numerically. For 
A = (2.3A)"/3o where A = 0.25 GeV such that Qpeak = 0.575 GeV, and integrating up to 
0.75 GeV, we obtain Iq = 0.536,0.441,0.389 and 0.355 for n = 2,4,6 and 8 respectively. 
From multiplicity considerations the preferred value is Iq ~ 0.5 according to the estimation 
from eqn. (^.S|), such that the values are reasonable, though the n = 6 and 8 cases are 
possibly too small. 



The cluster mass distribution corresponding to n = 2, 4 and 8 are shown in fig. 19 
At first sight, it seems that the n = 8 distribution matches the HERWIG default numbers 
best, hence implying that the low-energy tail of the as for n = 4, or 2, is too high. On 
the other hand, if the tail is not sufficiently high, as in the n = 8 case, there would not be 
sufficient contribution to Iq. Hence there seems to be a contradiction here. However, one 



reason why the Gaussian approximation in eqn. ( 5.13 ) can rectify the small-mass oriented 
skew of the cluster mass distribution is that there is some contribution from above 
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Figure 19: The cluster mass distribution generated using the as derived from power corrected 
beta function. 

(0.75 GeV)^. In HERWIG, the emission cut-ofF is implemented as being controlled by the 
gluon effective mass rather than as having a sharp cut-off, so that there is inevitably such 
a contribution. The proper study of this, taking account of the mass effects, is beyond 
our intention, and is presumably best left to a study incorporating experimental data, if 
possible using a generator that is based from the outset on the modified-as prescription 
rather than a modification of HERWIG as presented here. 

One way to study the g ^ qq dynamics might be to measure the compensation of the 
SU(2) light flavour charge. This should occur close to the electric charge compensation. 
Let us now turn to a discussion of the analyticity of low energy as- 
As stated above, the as defined by eqn. (^) has complex poles. This unwelcome in 
the sense that it forbids a spectral representation, viz: 

as{Q') = - r^^Psif^'), (6.8) 

Jo p^ + 

where psip-'^) is the spectral density function. As as is related to the gluon two-point func- 
tion, the physical interpretation would be that in the confined vacuum the spectral density 
associated with the colour-octet state is ill-defined. On the other hand, colour-singlet two- 
point functions should have a valid spectral representation such that the complex poles 
must somehow cancel. We add that in our convention, a physical spectral density function 
is negative definite. However, if we impose this condition, we see that as must continue to 
grow towards low Q^. 
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We may rephrase this statement as follows. If colour is confined in the large distance 
limit, an Q5 that describes the interaction of the unconfined gluon must vanish in the soft 
limit. But this is in violation of the spectral representation if the spectral density function 
is negative definite, such that there must be complex poles. These complex poles must 
cancel in quantities that are associated with colour singlet states. 

Let us denote the above as, eqn. ( |6.5| ), that describes gluons both in the confined and 
the asymptotically-free vacuum, by af"'- We then write the nonperturbative physics that 
cancels the poles by a^^' . Then for the overall process of hadronization, we have: 



h 

as 



»lu. 



I had. 

+ as 



(6.9) 



The former term is universal, but the latter term is not necessarily so. We also note 
that from the condition of finite emission probability, we required the first term to vanish 
as a power towards low . This is natural considering that an unconfined gluon can 
not propagate in the confined vacuum. On the other hand, as the concept of emission 
probability does not hold much meaning in the hadronic phase, there is no such requirement 
governing the latter term except that it is confined to the region of small Q^, and it would 
be reasonable for it to remain non-zero. 

Now let us shift contributions between the two terms such that there is no complex 
pole in either of the as, but the two terms, or more strictly their moments, are affected 
minimally, in some sense which we shall discuss later, for positive and real Q^. Then we 
can write: 



«5 



"5 +«5 ■ (6-10) 

The former 'confinement' term is still at least almost universal, a^'^^^ is expected to be 
the contribution to the observable under attention from the decay of the nonperturbative 
clusters and hadrons, and is expected to be effective at a lower scale than the former term. 

Using a^ which has a dispersive representa- 
tion, as do its two constituents a™"^^' and a^'^^^ , 
we can calculate power corrections following the 
procedure of ref. [18|. When doing so, the quanti- 



scheme 




F2,l 


-^3,1 


n = 2 


0.378 


0.394 


0.381 


n = 4 


0.394 


0.419 


0.401 


n = 6 


0.390 


0.421 


0.404 


n = 8 


0.385 


0.420 


0.404 


Webber's as 


0.511 


0.450 


0.410 



ties of main interest are the first few moments of 



Table 1: First three moments, Fp_i(2 
GeV), of low energy as- 



as, I.e., 
by: 



Fp,i{Q), where Fp^q{Q) |TH] are defined 



dk 



(6.11) 



a^g"" has complex poles and therefore can not be used to calculate power corrections 
directly in the dispersive approach. However, if our argument above is correct, we can 



a^g^' . As these moments can be calculated numerically for any choice of 
as so long as it is integrable, it then follows that one may calculate them for a^"' as 
an estimate of the moments of a^°°^'. We take the n = 2,4, 6 and 8 cases with A set to 
(2.3A)"'/3o and compare the result with as of ref. |15]. For Q = 2, the result is shown in 
tab. 0. 



choose a™°^' 
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The as of ref. [15| is defined to describe event shapes and is therefore an a^. Physically 
we would expect that both a™"^^' and a^^'^^^ would be positive definite. If so, we have 

> a™'^^'. A more important criterion is that this is satisfied in the moments. We see 
from tab. |^ that this is indeed the case. Furthermore, we see that a^'^^^ contribution only 
accounts for 0(20%) of the contribution to the moments. 

From the requirement that a^g'^^^ mainly affects the region of lower than a^°'^^-, or 
more strictly oi^g^' , we expect that the higher moments of and a™"^^' are similar. This is 
also satisfied in tab. |^. However, the similarity is exaggerated by the choice of the cut-off 
scale, 2 GeV, which is high compared with the hadronic scale. 

This argument governing the moments also implies that the moments Fp^i{Q) derived 
from experimental data must be greater than the contribution from q™"^^', such as those 
listed in tab. |l] except the one in the bottom row. From the data compiled in ref. pO| , 
which we reproduce in tab. ||, -Fi,i(2 GeV) is in the range 0.391 to 0.560. Although 0.391 
measured in the jet broadening observable is quite low, it is not in contradition of this 
statement. 

Let us proceed to estimating the ambiguity of Fi^i 
in the process of analytization in between eqns. (|6.9|) and 



( |6.10 ). We can do this by simply subtracting off the com- 



plex poles in a|'" and calculating the corresponding contri- 
bution to Fi i. This analytization procedure is equivalent 



to that adopted in ref. [21|. 

We first define for u and v real: 



Variable 


Fi,i(2 GeV) 


Bx 


0.4508 ± 0.0225 


Bw 


0.3911 ± 0.0305 


l-T 


0.4976 ± 0.0087 


C 


0.4527 ±0.0110 




0.5602 ± 0.0224 



u + iv 



A 



(6.12) 



Then Q / K is given by one of the roots: 



Table 2: Moment Fi,i(2 GeV) 
of low energy as, obtained by fit 
to event shape variables. Num- 
bers taken from ref. 120(1. 



^ M/(/3oA" 
A^ \ u + iv 



2/n 



a 



U + IV 



(6.13) 



where we defined for convenience a = A/{PqA'^) and = n/2. Let us specialize in the case 
of even n. It turns out that out of the u roots of eqn. ( |6.13| ), the only relevant one is the 
one whose phase is 1/z^ times the phase of n — iv, where this phase is taken to be between 
— vr and +11. 

At the pole we have: 



a. 



/?olog 



A2 



+ 



2A 



log 



a 



u + iv 



+ 



U + IV 



V 



(6.14) 



Hence: 



from which we derive: 



u + iv = a exp {u + iv) 



V 



tanv' 



e"- smv 



(6.15) 
(6.16) 
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Contour of poles for a§ on complex (A^qQ") plane 
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Figure 20: The contour of poles for on the complex A/(/3o(5") plane. 



As a is positive, we require that v/ sinv is also positive. With this condition, the contour 



of poles in the u + iv space is plotted in fig. 2C. a tends to +0 as — > +oo and +oo as 



u 



-oo. 



Eqns. ( 6.16| ) are consistent with eqn. ( |6.14 ) only if: 



\v\ < UTT. 



(6.17) 



Thus for n = 2,4, there are two solutions which are complex conjugates, whereas for 
n = 6, 8 there are four solutions, and so on. We note furthermore that the lowest solution 
is complex only if a > e~^. 

The poles of as can be found for given A and hence a by solving eqns. ( |6.16| ) to obtain 
V and u in terms of a. The residue on the plane is given in terms of the beta function 
evaluated at the pole, /3 = /3o(l — u — iv), as: 

(6.18) 

at pole 

When there is no complex pole, the spectral representation of eqn. ( |6.8D has the well- 
known solution in terms of the discontinuity of as on the negative real axis: 

Ps{^^') = 1^ [as{^^'en - «s(//'e— )] . (6.19) 



residue ^ 
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As stated earlier, a physical ps is negative definite in our convention. The generalization 
of eqns. ( |6.^ ) and ( |6.19| ) to include contribution from the complex poles yields: 



;U2 + Q2 



E 

pole 



(1 - Q2/^2 



(6.20) 



jj,^ at pole 



This general expression does not depend on the particular form of as but assumes that 
the complex singularities are all simple poles, as can be analytized by subtracting the 
complex poles, or by adding: 



pole 



(3{p^) (1 - Q2//z2^ 



(6.21) 



at pole 



Applying the formulae to the as in eqn. (|6.5| ), the resultant analytized as has a negative 
definite spectral density, given by: 



1 

To 



(6.22) 



n is assumed to be even. Denoting the analytized as by ZSg, as ots continues to grow 
with decreasing Q^, as{0) is a finite number, which can be calculated either by integrating 
eqn. ( |6.22 ), or more simply by summing eqn. ( |6.21| ) at = 0, i.e., summing 1//3 over the 
poles. 

'as{Q'^) has an unphysical 1/Q^ behaviour at large Q^, which must be cancelled by 
additional contribution from ps- This additional contribution can be defined as a negative 
definite quantity only if the sum of the residues Q'^/P, eqn. ( |6.1S| ), is positive. By evaluating 
eqn. (|6.18| ) in the A ^ oo limit, we see that this is always satisfied when u = n/2 is even, 
but violated for odd i' when A is large. 

The ambiguity of Fi^i is estimated by integrating eqn. ( |6.21 ) according to eqn. ( 6.11J ). 
We obtain 5Fi^i = 0.071,0.016 and —0.052 respectively for n = 2,4 and 8. A, as before, is 
given by (2.3A)"/3o. 

We see that the ambiguity, amounting to up to ~ 10% of Fi^i, is under control. Had 
we started from the unmodified perturbative as |21|, the corresponding SFi^i would be 
undefinable. 

As the procedure followed here gives an estimate only, it is not possible to draw defini- 
tive conclusions based on these numbers, but comparison with tabs. |l| and ^ shows that the 
magnitude of SFi^i is consistent with data, and in our opinion + 5Fi^i, corresponding 
to 05, also does not show disparity with tab. || that renders any particular choice of n 
completely unrealistic. 

We have so far calculated the confinement contribution of eqn. ( pTOD to and the 
ambiguity due to analytization in between eqns. (|6.9| ) and ( 6.10| ). It is not possible to 
calculate the decay contribution within our present framework. A very naive expectation 
would be that since the decay process occurs at lower energies than the confinement process, 
the highest energy component of a^'^'^' would be given by the same 6as that cancels the 
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complex pole in a^"'. If we accept this, we would expect that although is non- 

universal, its size is estimated by 5Fi^i that we have calculated above. Again, all cases are 
consistent with data and is in agreement with the small size of the a^^*^' contribution to 
Fi 1 seen in tab. ^ or in the difference between tabs. |^ and ^ 

One may ask how our results depend on the input parameters. We recall that we 
determined our value of A by tuning the peak position of as, as in eqn. (|6.7|). We may shift 
the peak position slightly and see how this affects our prediction. For n = 4, when Qpeak 
is modified from 0.575 GeV corresponding to Qmm = 0.4 GeV to 0.5 GeV corresponding 
to Qmin = 0.25 GeV, which we consider to be unreasonably small, 5-Fi,i decreases from 
0.015 to —0.004, while Fi^i from integrating a^'^' itself increases from 0.394 to 0.424. We 
have also calculated the case Qpeak = 0.65 GeV such that Qmm = 0.55 GeV which again 
is unreasonable. This yields 6Fi^i = 0.032 whereas Fi^i is reduced to 0.369. We observe 
that -Fi^i + 5Fi^i is less sensitive to A than Fi^i itself is, and in any case the dependence is 
small enough to trust the numbers calculated using the procedures presented here. 

7. Conclusions 

We have reconsidered the phenomenon of local charge compensation from the viewpoint of 
colour-preconfined models of hadronization. 

We have suggested that the study of minijet charge as a function of the resolution 
parameter ycut oc could provide information on the nature of dynamics governing 

hadronization. In particular, the scale at which local charge compensation is maximally 
violated, i.e., where the minijet charge is peaked, is the scale at which nonperturbative, 
confinement, dynamics sets in. We have demonstrated by simulations using HERWIG that 
modifications in the parameters governing hadronization could lead to significant differences 
in the behaviour of the minijet charge. We have made a comparison with PYTHIA and 
found that the predictions of the two generators are similar. 

The values for minijet charge depend on the scheme used for clustering minijets. Our 
default procedure uses the Cambridge algorithm, but we have demonstrated that the com- 
bination of the JADE algorithm in the confined region and angular-ordered Durham al- 
gorithm in the perturbative region results in less contamination from misidentified tracks. 
We have proposed a 'preclustered' scheme which shows marked improvement compared 
with the other algorithms when the preclustering scale 7/had is tuned to (0.5 GeV)^. 

In addition to the e^e~ case, we have presented a simple analysis of the minijet charge 
at hadron collisions and found that due to the extra contamination in this case, the peak in 
the minijet charge is weakened. The peak disappears completely when the soft underlying 
events are added according to the default option of HERWIG. 

We have studied a modification of the HERWIG default procedure to express the 
cluster-splitting dynamics at the end of the parton shower phase in terms of emissions 
that are considered unresolved in the course of the ordinary parton shower. Through this 
analysis, we have found that the lowest scale of cluster-splitting dynamics is about 0.4 
GeV. 
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This analysis also suggests that because the product of neighbouring cluster masses is 
smaller in this case compared with the HERWIG default procedure, the peak height in the 
minijet charge is necessarily smaller than is suggested by HERWIG. On the other hand, 
the relaxation of HERWIG's sudden transition to the confinement phase is expected to 
smear this peak. 

Our simulations for this part of the study were carried out using a simple modification 
of HERWIG with several approximations. It is hoped that simulation using Monte-Carlo 
event generators of the next generation, that are based from the outset on the modified- 
as prescription, will also be available in the future. We emphasize that comparison with 
experiment can nevertheless be made using our present approach to yield useful information 
that has less theoretical prejudice, but once a particular model of as is available, such as 
the one proposed in this paper, simulation based from the outset on modified-a^ could 
leave less ambiguity in the theoretical interpretation. 

We have discussed the possible implications of our findings to the underlying nonper- 
turbative dynamics. In particular, we studied a model in which the transition between 
the perturbative phase and the nonperturbative phase during the course of hadronization 
is driven by the semi-perturbative splitting g — > qq, possibly enhanced by the string-like 
confining force. This gives rise to an as which vanishes as a power towards low but 
has complex poles. We have argued, in agreement with available data, that the moments 
of this as gives an estimate of the part of the power corrections to event shapes that is 
universal and comes from soft gluon emission, whereas its complex poles give an estimate 
of the ambiguity, which is found to be small. In this picture, the remaining part of power 
corrections is related to the cascade decay of the nonperturbative, confined, clusters, and 
in our present framework can only be estimated to be of the same order as the ambiguity 
of the former part. 
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